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In the pedestal of a tokamak, the sharp radial gradients of density and temperature 
can give rise to poloidal variation in the density of impurities. At the same time, 
the flow of the impurity species is modified relative to the conventional neoclassical 
result. In this paper, these changes to the density and flow of a collisional impurity 
species are calculated for the case when the main ions are in the plateau regime. In 
this regime it is found that the impurity density can be higher at either the inboard 
or outboard side. This finding differs from earlier results for banana- or Pfirsch- 
Schliiter-regime main ions, in which case the impurity density is always higher at 
the inboard side in the absence of rotation. Finally, the modifications to the impurity 
flow are also given for the other regimes of main-ion collisionality. 

I. INTRODUCTION 

Conventional theory of neoclassical transport in tokamaks is not applicable to 

regions where the pressure and temperature profiles are very steep, such as the pedestal 
at the plasma edge. As the radial scale length decreases, poloidal variation arises in the 
temperature and density of each species. In an impure plasma, typically the first quantity 
to develop a poloidal variation is the impurity density, and indeed,strong poloidal impurity 
asymmetries have been observed in experiments js), Q|. In Refs. js - 7 1 , neoclassical theory for 
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an impure plasma was extended to allow for larger gradients than are usually considered. 
Specifically, the gradients were allowed to be so large that the friction between the bulk 
ions and heavy impurity ions could compete with the parallel impurity pressure gradient, as 
is typically the case in the tokamak edge. Mathematically, this means that the parameter 
A = bvaz 2 was assumed to be of order unity, but the poloidal Larmor radius of the bulk 
ions divided by the radial scale length associated with the density and temperature profiles 
5 = pe/L± was assumed to be small. Here z is the impurity charge number, = Ln/\ is a 
measure of the ion collisionality, \ is the bulk ion mean-free path, and L\\ is the connection 
length. It was shown that the impurity dynamics then become nonlinear, and if the pressure 
and temperature gradients of the main ion species are sufficiently steep, the impurities are 
pushed to the inboard side of the flux surface. 

Recently, the in-out density asymmetry A = nn/n-L was measured for boron impurities 
in Alcator C-Mod [3-]. Here, n# and refer respectively to the impurity density at 
the high-field- side midplane and low-field- side midplane of a given flux surface. It was 
observed that A could be either less than or greater than one. A comparison was made 

n 

to a theoretical model of impurity asymmetry in strong gradient regions [7J in which the 
primary ion species was assumed to be in the Pfirsch-Schliiter regime of collisionality. This 
model predicts that A must be more than one, and for the parameters of the Alcator C-Mod 
experiments, the predicted A was systematically closer to unity than the measured ratio. 
One factor which likely contributes to the discrepancy is that much of the data were taken 
in a region in which the main ions were in the plateau collisionality regime rather than 
the Pfirsch-Schliiter regime. Reference js} therefore suggests that an analogous theoretical 
model should be developed for the plateau regime, and it is the purpose of this paper to 
present such a model. Impurity asymmetry in the banana collisionality regime has been 
analyzed previously in 0, (J. Other than the collisionality, the present work uses the same 
orderings as the previous models: A ~ 1, z ^> 1, and 8 <C 1. 

The poloidal rearrangement of the impurities affects the impurity velocity due to the 
requirement of mass conservation. In previous work on the banana and Pfirsch-Schliiter 
regimes, this alteration to the impurity flow was not explicitly calculated. However, pedestal 
impurity flows are measured routinely in experiments 3|, |4( , so impurity flows represent an 
important point of comparison between experiment and theory. The measurements and 
conventional neoclassical theory often disagree. In particular, when the main ions are in the 
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plateau or banana collisionality regime, the measured impurity flow is greater in the direction 
of the electron diamagnetic velocity than predicted. Consequently, in this paper we give 
explicit forms for the modified impurity flows, and we examine whether the modifications 
are sufficient to reconcile neoclassical theory with the experimental measurements. 

The remainder of the paper is organized as follows. In Sec. [XT] we describe the kinetics 
of main ions in the plateau regime. In Sec. II III we analyze the parallel momentum equation 
for the impurities and derive an equation that governs their poloidal rearrangement. We 
show approximate solutions in several limits and numerical solutions are also presented. In 
Sec. HV] we explore the modification of the poloidal impurity rotation due to the presence 
of large gradients, discussing all regimes of main-ion collisionality. Finally, the results are 
summarized and discussed in Sec. |V] 



II. KINETICS OF MAIN IONS IN THE PLATEAU REGIME 

The plasma is assumed to consist of hydrogenic ions (i) in the plateau regime, collisional 
(Pfirsch-Schluter) impurities (z), and electrons (e). The calculation does not depend on the 
collisionality regime of the electrons. The magnetic field is represented as B = /(■0)V^ + 
V93 x Vip, where ip is the toroidal angle and 2irip is the poloidal flux. Throughout this analysis 
we will use a poloidal angle coordinate d which is chosen so that B ■ V?9 is a flux function. 
This coordinate makes flux surface averages convenient to evaluate ((Y) = (27r) _1 J^Yd'd 
for any quantity Y), and this coordinate is equivalent to the $ used in 5K7|. We assume a 
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model field magnitude b 2 = 1 — 2ecos?9 where b = B/ (B 2 ) and e = r/R is the inverse 
aspect ratio. We must assume e C 1 from the beginning of the analysis in order for a plateau 
regime to exist. 

The gyroaveraged ion distribution function in the plateau regime is then given by [h| 
ft = fMi + In where 

is a stationary Maxwellian and a flux function, 

fn - -fMi-jr + H t - f Mi v\\^ y- + ^r + ^r-)i ( 2 ) 
Pi = n i0 Ti, fij = eB/rrii is the ion cyclotron frequency, primes denote d/dip, $0 = 
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$ x = $ - $ 0) 

i>i sin i? — X|| cost? ^ I" , . „ cos$ 
2~m + f,- 



suit? 



z>j = viqR/vi is the normalized collisionality, x = v/vi, Vi = (2Ti/m i ) 1 ^ 2 



(3) 



Qi = fMi 6 -^ [(2x1 + xi)(2x 2 - 5) + yb\2x\ - x\)] , (4) 

and y is a velocity-independent coefficient which will be determined by the requirement 
of ambipolarity. In a pure plasma this requirement leads to y = 1, but the presence of 
impurities will alter the value. Also, y must be a flux function so that V • (n^Vi) = 0. 

III. IMPURITY DYNAMICS 

The parallel momentum equation for the impurities is taken to be 

= -zn 2 eV||$ - T i V\\n z + R z q (5) 

where R z m is the impurity-ion friction. The parallel viscosity of the impurities has been 
neglected since it was shown in Ref s| to be smaller than the pressure gradient if 5/zi>a <C 1, 
which is usually the case in the tokamak edge. As also shown in that paper, the impurity 
temperature is then equilibrated with the bulk ion temperature and is therefore constant 
over the flux surface. The poloidal electric field — V||$ can be obtained from the quasi- 
neutrality condition zn z = n e —rii using n e — (1 + e<fe\/T e )n e o(ip) and using the distribution 
function (jSJ) to calculate the ion density: 

Hi = n i0 ( 1 - — + eN s sin ■& ) , (6) 



where 



The result is 



Ti 2T-n 

where 2n /T = n e0 /T e + n^jTi. Equation (jSJ) then becomes 



V|| (z 2 n z + zn i0 eN s sin 0) , (8) 



l + an)V\\n+—— Vy [N s smti ) = —— 9 
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where n = n z /{n z ) is the normalized impurity density and a = (n z ) z 2 T j (2noTj) . In the 
rest of the analysis we will order a ~ 1, which is equivalent (for T e ~ T 4 ) to the ordering 
z e// -1-1. 

Next, the ion-impurity collision operator Cj 2 is inserted in = — m,j J d 3 vv\\Ci z to 
write 

= -J d 3 vm lVll u tz (c(fi - U) + ^^ii/ioj (10) 

where 

„ 2v \\ 9 , <9 

is the Lorentz pitch-angle scattering operator, A = w^/(i?w 2 ), = 37r 1 / 2 /(4rj 2 a; 3 ), and 
r iz = 3(27rTj) 3//2 egm| //2 /(n 2 ^ 2 e 4 ) is the ion-impurity collision time. To ensure V- (n z V z ) = 0, 
the parallel impurity flow velocity must have the form 

B n z 

where K z (ip) is proportional to the poloidal velocity. Using the main- ion distribution func- 
tion ([2]) we then obtain 

d 1 ( , , yb 2 n i0 T(\ mini K z B 

R -^-^{ p - + ^)'^r +Q - (13) 



where 

Q r = nrii / d 3 vu iz V\\Qi 



Xw 

tt6(x\\) siwd k — - — k cosi!) 

11 x 2 + Pf 



(14) 



For z>j — > the integration results in 



Qr = 3 ^oIT[ cosr& _ (15) 

To rewrite Eq. (Q in dimensionless form we introduce the ratio of the temperature and 
pressure scale lengths r) = piT- / '(TipQ , 

m i l v'i (IQ) 

and 

_ ^/TTzT T iz n z ViB ■ W 

where B = (B 2 ) 1 ^ 2 . Notice that g, r*, and (T iz n z ) are "^-independent, and the formal 
magnitude of r* ~ (zz/,,)" 1 has not yet been fixed. Equation fl2]) now becomes 

(l + am)^ = ^|n+^-07[3 + (l + y)rJnoosi? + iC ( ^^} (18) 
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where we have used d(N s sm , d)/d'd ~ (N s ) cos$. (Other terms of order e have already been 
discarded in deriving the distribution function fl2]).) Integrating Eq ( Tl8]) over $ yields a 
solubility constraint which can be used to determine the poloidal impurity rotation, 

K z = {-1 - f (nb 2 ) + [3 + (1 + y)n] e V (n costf) } , (19) 

and Eq ( Fl8|) becomes 

n _ b 2 + m^^ n _ ( nb 2^ 

+ [3 + (1 + y)n] er] (b 2 (ncosti) -ncosti)] . (20) 



The cos $ terms above can be significant despite being proportional to e, for the other drive 
in the equation is the variation in b, which is also 0(e). 

To make further progress we will calculate the coefficient y by requiring ambipolarity. 
Due to the smallness of the electron mass, the ambipolarity condition is approximately 
Tj = —zY z . As in the conventional plateau-regime calculation for a pure plasma, the main- 
ion flux is 

r, s (r , . v,) = ^g^'^ fa - 1) CD 

where Qio = eBo/rrii. The impurity flux is driven by the impurity-ion parallel friction force 

lRzi\\ 



(T- ■ W) = < ) , (22) 



where R Z i\\ is given by Eq. (I13p with K z from (I19p . We find 

+er] 1 3 - [3 + (1 + y)n] (n cos i9) }) . (23) 

The condition for ambipolarity then gives 

z<?T*or x +r]- 1 ((n/b 2 ) - 1) - 3e(n6~ 2 cos$) + (3 + r*)e (ncos$) 
y= ze 2 ^^ 1 + 2-!((^ 2 ) - 1) - en (ncostf) ' ^ ' 

The pure plasma limit y — 1 is recovered as a — > 0. 

The system (1201) and ( 1241) describes the poloidal rearrangement of the impurities. While 
(1201) is similar to the equations found if the main ions are in the banana 



6j or Pfirsch- 

Schluter regimes [3] , (120]) has several different terms, and also the radial scale length entering 
g is different (i.e. only the pressure scale length appears, rather than a combination of the 



pressure and temperature scale lengths.) As in Refs.js-7], g measures the steepness of the 
bulk ion pressure profile. In conventional neoclassical theory g is assumed to be small, 
which implies that the friction force is smaller than the parallel pressure gradient. We next 
examine how the integro-differential equation ( 1201) can be solved analytically in a number 
of limits. 

a. Weak density variation. If n — 1 ~ 0(e) then we can expand n = 1 + n c cos$ + 

n s smi9 + 0(e 2 ) with n s and n c both ~ O(e). The solution of Eq (120]) is then found to be 

2-77[3 + (l + y)rJ 

n s = eg(l + a)— — 25 

(1 + a) 2 + 3^(1 + riy/2) 2 

n c = -eg\l + V y 2 — [ ^ g; * J 26 

It can be noted from these expressions that as p£ becomes larger, the impurities first develop 
an up-down asymmetry and then an in-out asymmetry. This same behaviour is found in 
the banana and Pfirsch-Schliiter regimes. However, in the plateau regime the asymmetry is 
proportional to the new factor 2— r)[3+(l+y)T*\, which means that the sign of the asymmetry 
can be changed depending on the magnitude of rj, r* and y. If rj > 2/[3 + (1 + y)r*], the 
impurities will be pushed to the outside of the flux surface. This result is different from the 
analogous n— 1 ~ 0(e) 1 limits when the main ions are in the banana or Pfirsch-Schliiter 
regimes. In these cases, in the absence of rotation, the impurities were pushed to the inside, 
regardless of the ratio of the pressure and temperature gradients. 

b. Large gradients. In the g > 1 limit, corresponding to a large pressure gradient, we 
can expand (1201) in g^ 1 . To lowest order, the right-hand side of ( 120]) must vanish, giving 
n ~ nl (h) where 

b 2 

n = (27) 

1 + (r]y /2)b 2 -er][3 + (l + y)n] cos ^' V ; 

In this case there is only in-out asymmetry. Expanding in e then gives 

n = l + 2e(,S-l)cos^ (28) 

where S — T)[3 + y + (l + y)r«] / (2 + rjy). (This same result can also be obtained by a g S> 1 
expansion of (126 p .) The impurity density evidently may be higher at either the inboard side 
(S < 1) or outboard side (S > 1). This finding too differs from the corresponding 3 > 1 
limits when the main ions are in the banana or Pfirsch-Schliiter regime. In these cases, 
the impurity density is always greater at the inboard side (even when there is significant 
rotation). 
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c. Numerical solution. For a « 1, equation (1201) may be solved numerically with 
the following iterative procedure. A small number (5-10) of poloidal Fourier modes are 
considered. An initial guess for is used to compute y and the nonlinear term an dn/dd. 
An improved n{d) is then calculated using ( 1201) . and the process is repeated until convergence 
is achieved. Typical results are shown in Figure [TJ Figure [2] shows the in-out asymmetry 
factor 

n{$ = 7r) 



.4 



(29) 



n{d = 0) 
over a wide range of parameters. 

Figure ([3]) shows the in-out asymmetry A for e = 0.3, g = 10, and the trace impurity 
limit y — > 1. A nearly identical plot can be generated using the g ^> 1 expressions ( 1271) or 
( 12 8 j) . although the precise value of A in the A < 1 region is somewhat different due to the 
fact that e = 0.3 is not much smaller than one. 

ri=0.2 (b) ri=0.7 



(a) 





0.5 



— g=10 

g = 3 
g=l 



FIG. 1: (Color online) Normalized impurity density as function of poloidal angle, calculated by 
numerical solution of ([20]) and (p4^) . The parameters used are e = 0.3, r* = 0.5, z = 5, and 
a = 0.25. Other values of a from to 1 produce nearly indistinguishable results. 



IV. POLOIDAL IMPURITY ROTATION 

If the impurity density varies on a flux surface, the impurity poloidal rotation will be 
different from the one derived in conventional neoclassical theory. Using (fT2|) and (fT9|) . we 
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(a) t* = o.i (b) x. = i 




FIG. 2: (Color online) Contours of the in-out asymmetry A, decreasing monotonically with r/, for 
(a) r* = 0.1 and (b) r* = 1.0. The other parameters are e = 0.3, z = 5, and a = 0.25. Results 
for a = are nearly indistinguishable. Solid contours run from A = 3 to 1 in steps of 0.5. Dashed 
contours decrease from A = 0.9 to 0.6 (a) or 0.2 (b) in steps of 0.1. 




FIG. 3: (Color online) Contours of the in-out asymmetry factor A in the g> 1 and trace impurity 
(y — > 1) limit and with e = 0.3. Solid contours range from 3.5 to 1 with spacing of 0.5, and dashed 
contours range from 0.9 to 0.1 with spacing of 0.1. 
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can write 



where 



X 



ypl 



1 + 1 



B»K Z 



-X 



IB, 



n 



ne {B 2 ) 



Tj dn i0 3 dTj 
n i0 dip 2 dip 



1 + y (nb 2 ) - [3 + (1 + y)n) erj (ncostf) 



(30) 



(31) 



is constant on a flux surface. The definition of X was chosen above so that in the trace 
impurity limit (a — > 0, y — > 1) and if n z is also uniform on a flux surface (i.e. g — > 0), then 
X — > 1. This limit reproduces the conventional neoclassical result 

an 

Figure H] shows the scale factor X for various values of n, r#, and g. The figure was 
calculated using e = 0.3 and a — > 0. It is evident that when </ > 1, the poloidal flow 
can be significantly suppressed compared to the conventional neoclassical result if rj and 
approach one. The situation is only slightly different when the relative impurity strength a 
is nonzero, as shown in Figure [U This figure is equivalent to Figure Hla but with a raised 
to 0.25 and z = 5. When r* C 1, the flow now becomes slightly enhanced compared to the 
conventional neoclassical result. 

When the main ions are in the banana regime, the poloidal impurity flow can be calculated 
using the K z derived in Refs. fl and Q. The result is 



V. 



ban 



B, 



-u 



1 



ban 



(32) 



where 7 " an = eL ±Mn (B 2 ) u/T t and L^ ban 
impurities and large aspect ratio, 



I {Pi/Pi - (3/2)7y/T;). In the limit of trace 



-0.33/ c 



and 



fc 



3(B 2 ) 
4 



Ac 



/ dTi 
e{B 2 )lkp' 

XdX 



n 



VI - XB) 



is the effective fraction of circulating particles. Therefore, in this limit, 



V. 



ban 

zti 



IB, 



ne (B 2 ) 



Ti dn 



t0 



+ (-- + 0.33 /c <^ 2 ))^ 



dTi 



(33) 



(34) 



(35) 



n i0 dip 

The expression for u in various other limits (arbitrary aspect ratio and high level of impu- 
rities) is more complicated and is given in Ref. js]. 

When the impurity density is nearly constant on a flux surface, (1541) gives the conventional 
result f c ~ 1 — lA6y/e. For insight into how f c is modified when the impurity density varies 
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FIG. 4: (Color online) The factor X which scales the poloidal impurity flow in the plateau regime, 
calculated for a = 0. The horizontal axis is the same for all plots. Contour spacing is 0.05. 
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significantly on a flux surface, consider the limit n = <5 (7? — n) in which the impurities are 
strongly peaked on the inboard midplane. Then (ny/\ — A£>) = \/l — \B max so f c ~ 1 — 2e. 

Similarly, when the main ions are in the Pfirsch-Schliiter regime, the poloidal impurity 
flow V^ s can be calculated using the K z derived in equation (26) of Ref. Q. For trace 
impurities, V^ s is found to be 



ne (B 2 ) 



Ti dn i0 _ _ , dTi 



2.8 <n6 2 )^ 
n i0 dip dtp 



(36) 



It was found in Refs. Q and Q that when the main ions are in the banana or Pfirsch- 
Schliiter regimes, the impurities tend to accumulate on the high field side, so (nb 2 ) > 1. In 
both regimes, this change decreases the signed V z $, shifting the poloidal impurity flow in 
the direction of the electron diamagnetic velocity relative to the conventional neoclassical 
prediction. We can model the impurity density variation as n = 1 — (A — 1)(A + 1) _1 cos^, 
implying (nb 2 ) = 1 + e(A — 1)/(A+ 1). As A increases above one, (nb 2 ) increases from one 
to 1 + e. For the banana regime, this increase in (nb 2 ) and the aforementioned increase in 
f c both lead to a decrease in the signed V z #, with the 0(y/e) increase in f c being the larger 
of the two effects. 

Note that in the method used in this section, the impurity pressure gradient p' z does not 
appear in the formulae for the poloidal impurity flow for any collisionality regime (as it does 
in, for example, equation (15) of 11]). In the conventional neoclassical formulae, the p' z term 
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is proportional to 1/z, so the term is formally small in our ordering. The absence of the p' z 
term is related to the fact that the impurity diamagnetic flow was dropped in Eq. f lT2|) in 
order to make the analysis tractable. 



V. CONCLUSIONS AND DISCUSSION 



In this paper we have investigated the poloidal rearrangement of impurities in the presence 
of large gradients for the case of background ions in the plateau collisionality regime. The 
calculation shows that when the temperature scale length is large compared to the density 
scale length (such that rj < 0.4 — 0.6), the impurities accumulate on the inboard side, whereas 
they accumulate on the outboard side in the opposite case. In standard tokamak operating 
regimes, t] < 0.5, so impurity accumulation at the inboard side is more likely. However, t] can 
be larger than 0.5 in the I-mode regime of Alcator C-Mod jljj, so the strong ^-dependence of 
A predicted by the theory may be experimentally testable. (Impurity asymmetry in I-mode 
has not been measured as of this writing.) 

One way in which the present calculation could be extended would be to account for 
the large radial electric field E r which arises in the pedestal. It is found experimentally 
that the radial electric field in the pedestal can be lar ge e nough to make the E x B drift 
comparable to (Bg/B)vi, and it was recently shown in 10j that under these conditions, the 
plateau- regime ion distribution function can deviate from Eqs. (JSJHj). Although it would be 
desirable to include this effect in the present calculation of impurity asymmetry, doing so is 
not straightforward, for the following reason. Terms in the ion distribution function of order 
(pe/ajefMi affect the impurity asymmetry calculation to leading order, as demonstrated 
by the term with a factor of 3 in Eq. (1201. which arises due to the Q r term in Eq. ([TBI . 
However, the ion distribution function in [l0] is only determined to order (pe/a)e /mi, and 
to consistently determine all 0(e) corrections, the ion distribution function would need to 
be found using the full linearized Fokker-Planck collision operator rather than a Krook or 
pitch-angle scattering model operator. 

In all regimes of main-ion collisionality, the poloidal rearrangement of impurities results 
in changes to the the poloidal impurity flow. These modifications to the flow are of interest 
because when the main ion collisionality is in the plateau or banana regimes in Alcator C- 
Mod, impurity velocity in the pedestal is measured to be greater in the electron diamagnetic 
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direction than conventional neoclassical theory predicts When the ions are in the 

plateau regime, the calculation in this paper shows the impurity flow should be multiplied 
by the factor X relative to the conventional neoclassical prediction (in which the flow is 
always in the electron diamagnetic direction.) To explain the observed flows, then, X must 
be > 1, which can occur for r t < 1 (as in Figure El) When the ions are in the Pfirsch- 
Schluter regime, we find the poloidal impurity flow is indeed increased in the direction of 
the electron diamagnetic velocity due to the increase in (nb 2 ) above one. For banana-regime 
ions, the flow is shifted in the same direction due to both the increase in (nb 2 ) and also 
due to the increase in f c . However, the shift in the flow is also proportional to the small 
numerical factor 0.33 in (135 j) . so this effect is likely insufficient to explain the observed 
discrepancy between the measured and predicted flows. A different calculation, including 

;he E r effects discussed above but neglecting the impurity asymmetry, is discussed in Ref. 

4(; this calculation can also explain some but not all of the discrepancy. In future work, 
it may be possible to consistently account for both the E r and impurity asymmetry effects 
simultaneously to achieve better agreement between the calculated and observed flows. 
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